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1 Introduction 

The energy scale of cosmic inflation is largely unknown. We know that it must be higher 
than about an MeV from nucleosynethesis and less than the Planck energy scale. Various 
assumptions can tighten this range, but experimental handles are clearly needed. 

Currently, the most direct handle on the inflationary energy scale comes from (non)- 
detection of tensor fluctuations. This is because in typical models of inflation, the ampli- 
tude of tensor fluctuations is in one-to-one correspondence with the inflationary Hubble 
parameter, whereas scalar fluctuations are degenerate in H and the slow-roll parameter e. 1 

However, this conclusion has recently been challenged [2-5]. These authors have 
pointed out that although the amplitude of vacuum tensor fluctuations uniquely probes 
the Hubble scale, it is conceivable that dynamics involving the metric could produce large 
correlations. Concretely, these authors have proposed that a spectator x held (or fields) 
rapidly excited out of vacuum during inflation could produce significant radiative correc- 
tions to the tensor power spectrum. 

The purpose of this paper is to clarify precisely in what way such dynamics can affect 
the tensor spectrum. Previous work has consisted of order-of-magnitude estimates [3] or 
computations based on a perturbative expansion that differs from ours [2] ; we elaborate on 

1 This is typically parametrized in terms of the ratio r of tensor fluctuations to scalar fluctuations, 
the latter of which have been observed. This ratio is related to the scale of inflation by r ~ fOO x 
Vi n flation/(10 16 Gel/) 4 ; to date the non-detection of tensor fluctuations gives an upper bound r < 0.2, 
yielding V^tton £ 10" 2 M pi . [f] 
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these differences in the text. Here, we show how to organize both the production of x excita- 
tions and its effect on the tensor power spectrum in a systematic field theory computation. 
We develop an analytical method to compute the various perturbative contributions. 

We find that any effects from rapidly excited scalar fields on the tensor power spectrum 
is small, as long as perturbation theory holds. The reason is that in order for the field to 
become excited, its mass must change rapidly on a timescale At <C H~ l . Consistency of 
perturbation theory puts a lower bound on this timescale, restricting the non-adiabatic 
nature of the process. This leads to an upper bound on corrections to the tensor spectrum, 
rendering the corrections unobservable. 

The rest of the paper is organized as follows. In section two, we review the inflationary 
background dynamics and production of x- In section three, we introduce our formalism for 
computing the tensor spectrum and state the results; the detailed computation is relegated 
to the appendix. In section four we make some comments regarding an alternative scenario 
in which x has a nontrivial classical background, again concluding that significant effects 
on the tensors are generally ruled out. In section five we conclude. 



2 Background equations and x production 

Consider gravity minimally coupled to two scalar fields, 4> and %• The action is 

S = \ J ^x^g [M 2 R + d^d^ + d^X ~ 2V(& x)] • (2-1) 

We assume that the potential is such that the inflaton <j) has a large homogeneous part 
(j)(x, t) = Tp{t) + 5<p(x, t) which sources a period of accelerated expansion, while \ ac ^ s 
only as a perturbation. The gravitational background sourced by Tp is <?qo(*) = 1 = 
—a(t) 2 5ij, i.e. we assume a flat FRW metric. As usual, a(t) satisfies the Friedmann equation 
H 2 M 2 = V(lp)/3 and we take the de Sitter approximation a(t) ~ e Ht ; here and throughout 
M denotes the Planck mass. 

We want to study the effects of excited quantum fields x on the power spectrum of 
tensor perturbations. Inflation redshifts any such excitations away, so the initial conditions 
for field perturbations are the closest thing to a vacuum state that may be defined, the 
"adiabatic vacuum," made precise shortly. Each mode 

X(x, t) = j ^ e *-** k (t), xUt) = Xk(t)a k + X^)aL k . (2.2) 
is a damped harmonic oscillator with time-dependent frequency: 

Xk + SH Xk + u 2 Xk = 0, Co 2 (k,t)= (-^-) +m 2 (t). (2.3) 



a(t) 

Here we used rotational invariance to write the modes in terms of k = |k|. This equation 
can be mapped into an undamped oscillator by taking x ^ a 3//2 %, except that this oscillator 
will have frequency 

J = - V. (2.4) 
4 
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The idea proposed in [2, 3] is to give x a time-dependent mass such that it becomes massless 
at some point during inflation. At low k the mode's frequency will change non-adiabatically 
|o;/o; 2 | 3> 1 and the mode will be excited. We assume that m(to) = at some time to and 
take this to be a local minimum so that m 2 (t) ~ Th^t — to) 2 = (/c*/ao) 4 (t — to) 2 near to. 
Here clq = a(to) is the scale factor, and we have defined the (comoving) momentum scale 
which we will use to phrase all of our arguments. 2 
For any comoving mode k, the time-variation of lv is very slow away from to (fig. 1). 
Thus the adiabatic mode function 



xl(t) 



1 



y/2(2ir) 3 a 3 (t)u(k,t) 



exp 



|-*y dt'uj(k,t')^j (2.5) 



and its complex conjugate provide, to very good approximation, a basis of solutions to (2.3), 
normalized to satisfy the Wronskian condition, see eg. [9]. Before the production event, 
t <C to, we need to define the state of the x field. High momentum modes with k/a^> m 
should match the Bunch-Davies behavior x ~ 1/k. Non-relativistic modes, k/a -C m, will 
see the flat space vacuum, x ~ 1/m. The adiabatic mode function (2.5) reduces to both 
of these limits. According to the usual argument (the equivalence principle) we therefore 
take as an initial condition Xk = Xk an d set the state to be the one defined by au | 0) x = 0. 
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Figure 1. Left: Non-adiabaticity parameter |w/w 2 | as a function of time for various wavelengths. 
Right: Corresponding numerical time evolution of the x correlation function (xk(t)X-k(t)) = 
\Xk(t)\ 2 . Long wavelengths (blue) have |w/cj 2 | ^S> 1 near to (in this picture, to ss 300); these 
modes undergo a sharp non-adiabatic evolution, leading to interference between the positive and 
negative frequency mode functions, cf. (2.6). Shorter- wavelength modes have weaker effects (red, 
{; k i;,) or evolve completely adiabatically (yellow, k 3> &*). Notice that all mode functions have a 
nontrivial evolution beyond the usual redshifting: they scale like a; -1 / 2 , which is enhanced near t - 

From eq. (2.4), we see that the behavior of the modes of x near t = to depends on the 
wavelength. Very high momentum modes have real and slowly-varying frequency for all 

2 The prototypical model [2] achieves this up by a coupling g 2 X 2 (4>~ <t>a) 2 between \ an d the inflaton. In 
these models, k, — ao \J gp where <^ denotes the background inflaton velocity at t = to. This model suffers 
from large corrections to the scalar power spectrum, so here we just leave the mass of x as an unspecified 
function of time without positing its origin as a VEV of some other field. 
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times and thus evolve adiabatically. The more interesting case is that of modes with long 
wavelengths, k/aoH <C oo, so that the frequency evolves rapidly |w/w 2 | S> 1. These modes 
evolve non-adiabatically near to, (fig. 1). The timescale of this process is At = ao/fe*; we 
assume that this is shorter than the Hubble timescale, i.e. At <C H~ l . 

Sufficiently long after to, every mode again evolves adiabatically and so the general 
solution to (2.3) must be a superposition of the adiabatic solution and its conjugate, 



Xk (t) = a(k, t)x° k (t) + f3(k, t)x° k (t). (2.6) 

The coefficients a and (3 are the usual Bogoliubov coefficients; the Wronskian condition 
reduces to |a| 2 — 1/3| 2 = 1. The initial condition described above corresponds to a = 1,0 = 
as t — > — oo. Modes that evolve non-adiabatically pick up a non-zero (3 around t = to and 
then maintain this value for t > to- 

A late-time observer, viewing physics with respect to her naturally defined adiabatic 
vacuum, 3 would see the \ field modes excited with occupation numbers = \(3(k)\ 2 . In 
order to quantify the excitations of x> we thus want to find the value of j3 at late times. 
The simplest way is to solve the mode equation numerically, fig. 1. One can also find 
approximate analytic solutions and match these to (2.6), see eg. [2, 6]. 

Another approach, developed in [11], is to consider the analytical continuation of u to 
complex times. Later, we will generalize this approach to compute our loop diagrams. It 
can be shown that, if |w/w 2 | <C 1 for all time, then at any time t one can compute j3 by 
the integral [11] 

- L di IMh^{- 2i L dfu[k ' f) \ (2 - 7) 

Note that for |t — to| 3> At, ui/u is small and the large u makes the phase incoherent, so 
(3 is constant. The point is that near to this integral would naively pick up a non-trivial 
contribution; the problem is that the integral does not make sense in this region because 
the adiabaticity condition |w/u; 2 | <C 1 is violated. The solution is to adapt a trick from 
ID Schrodinger scattering problems: we can deform the contour of integration here off of 
the real axis in such a way that |w/w 2 | <C 1 along the entire contour. The frequency oj is 
a square root and thus has a branch cut on the complex plane, and as we go around this 
cut, the frequency picks up an imaginary part, so that the exponential picks up a real part. 
We give a detailed account of this method in the appendix. For the time being, the main 
point is that both the matching and analytical continuation methods have been used to 
find f3(t — > +oo), and they agree, yielding 



f3(t -> oo) = -exp <^ -- 



k \ 2 /3 aoH \ 2 



k* J \ 2 k* J 



(2i 



.e. the state defined by bk | 0} = where bk is the Bogoliubov-transformed annihilation operator. [10] 
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This formula is the key to our result that the effect on the tensor spectrum is small. The 
tensors are derivatively coupled to x but x is produced in the infrared, so loop effects are 
localized to a small and infrared portion of phase space. 

3 Computation of the tensor spectrum 

The basic observable quantity we are interested in is the two-point correlation function of 
the tensor fluctuations, defined by g /JiU (x,t) = g^t) + h fiu (x,t). Here g is the background 
inflationary metric described earlier. One typically defines the tensor power spectrum as 

A?(A;) = ^Pt(k) (3.1) 
where Pt(k) is the two-point correlation function of the tensor fluctuations, 5 

(2nfP t (k)5(k + k') = (3.3) 

s 

This is an in-in computation which we perform largely following [9]. We work in the 
interaction picture, so that all field operators evolve according to their free 6 equations of 
motion while the state evolves according to the time evolution determined by the interaction 
Hamiltonian. In other words, for any operator O, we have 

| O | ¥>(f) = (* | U-\u,t)0{t)U{U,t) | (3.4) 
where the explicit time-dependence of the operator is that of the interaction picture, and 

U(ti,t) = Texp^-i jf dt' H int {t')^ (3.5) 

where here Hi n t has both explicit time-dependence as well as time-dependence from writing 
all field operators in the interaction picture. The state | \P) is determined at the initial 

4 Previous studies neglected the expansion of the universe by simply setting H = 0. For subhorizon modes 
k 2> aoH this approximation is good as can be seen easily from the above formula since by assumption of 
non-adiabaticity At -1 = k* S> H; the derivative couplings of x to the gravitons mean that this is also a 
good approximation for our problem. 

J The Fourier conventions here are 

h ab (x,t)=J2j -^e ikx e S ab(k)h s k (t), h a k (t) = h B k (t)b s k + ht^)b s \. (3.2) 
We will work in the transverse-traceless gauge for the gravitons. The polarization tensors satisfy[12] 

£ab(k)e a cd ( — k) = S ac 5bd + S a dSbc ~ S a bS c d + S a bk c kd + 8cd.ka.kb 

8ackbkd 8 a dkt,k c 8b c k a kd 8bdkak c -t- k a kbk c k d . 

6 Here we are treating the x field as having a time-dependent mass; of course if this mass comes from eg. 
coupling x to Tp we still include this coupling in the "free" equations of motion. 
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time ij. In practice we are going to take tj — > — oo and set the state to be | 0), the state 
annihilated by the annihilation operators of all the relevant fields, as described above. 

In our problem the relevant part of the interaction Hamiltonian comes from expanding 
the Lagrangian for x minimally coupled to gravity, giving the interactions (see eg. [13]) 



C 



int 



-g 



(3.6) 



Here we have dropped the terms that vanish for transverse, traceless gravitons. The in- 
teraction picture fields have time-evolution governed by their free equations of motion. 
Ignoring the Planck-suppressed effects of x on the tensors, each polarization is a massless 
minimally coupled scalar field in de Sitter, 



Kit) 



H 



MV2k 3 



1 + i 



k 



aH 



(3.7) 



Meanwhile x has a more complicated evolution due to its interaction with the background 
inflaton field, as described above. The object of interest is (^|/i|/t?,|\I/)(i). From here out 
we drop the explicit reference to the state to maintain some sanity in the notation. Using 
(3.4), (3.5) and Wick's theorem, 7 we can expand this in a kind of Dyson series, in powers of 
the interaction Hamiltonian. Following the logic in [9] , it is possible to develop a Feynman 
diagram-style calculus to organize this computation. Unlike Weinberg we work in Fourier 
space directly, but otherwise the procedure is identical. 



k ■ s , 1 





Figure 2. Diagrams contributing to (hkhk) at leading order (1/M 2 ). The scalar lines are % 
correlators and the curly lines are gravitons. Later we denote the first diagram by A and the second 
by B. 

The standard result for the tensor spectrum comes from assuming that the metric 
fluctuations are in the Bunch-Davies state defined above and ignores interactions. In the 
language we are using here, this means it is just the tree- level correlation function: 




2M 2 k 3 



(3.8) 



At late times t — )■ oo this gives the usual scale-invariant spectrum of inflation. Corrections 
to this come from insertions of the interaction Hamiltonian in (3.4). In our problem this 
means that the first non- vanishing corrections will come at order 1/M 2 . 



7 Wick's theorem for time-ordered vacuum correlators is of course (TY\ i Oi(ti)} = ^2ij{TOi(ti)Oj(tj)}. 
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There are two diagram topologies contributing at first order, see fig. 2. Adapting 
the rules from [9] or just brute-force expanding (3.4), these can be evaluated in terms 
of integrals over time and loop momenta. Each vertex is labelled by a time argument 
which is integrated over from t — > — oo until the time of measurement t (in practice we 
measure correlators well outside the horizon, t — > oo; we have checked numerically that the 
loop corrections become very nearly constant after horizon exit of the "external" graviton 
momentum k, as with standard loop corrections in inflation). For each vertex, we may 
have an insertion on either the right or the left in (3.4), and we must sum over all of 
these possibilities. The loops of x fluctuations contain the information about x excitations 
because the propagators are enhanced as explained in the previous section. 

The first diagram, call it A, is a tadpole and is straightforward to evaluate. This dia- 
gram has been missed by previous treatments of the problem. Summing over polarizations 
s and adding the diagram with a right insertion to the one with a left insertion, one obtains 



s,R,L 



-i«5(k + k') 



dt'd 3 q 



f 4 (q,t')lm[hUt)hl(t')Y \ Xq (t')\ 2 



(3.9) 



where fi comes from the sum over polarizations. The second diagram B has four combi- 
nations of right/left choices. Summing over polarizations again, we find for the right-right 
diagram 



(hih s k/ ) RR = -5(k + k') / dt' 



dt" 



d 3 c 



(2tt) 2 



(3.10) 



x/ 3 (g,* / ,0^(*)^(<')^(*)^(Ox g (t')Xg(t")x,'(*')x,'(t") 



where k + q + q' = and / 3 is again from the polarization sum. The left-left diagram is 
just the conjugate of this, {h^Ji^/jLL = {hfJi^RR- ^he mixed-vertex diagrams contribute 
equally, 



Y,ChiK,) RL = +5(k + k') f df f dt" f 



d 3 q 
(2^f 



(3.11) 



x/ 3 (g/ > 0^(*')^(<)^W^(< /, )x g («')x 9 (<' / )x g '(*')x,'(<' / ). 
and (h s k h s k ,)LR = {h s k h s k ,) R L- Explicitly, the functions /3 and fi are given by 

/ 3 (<Z, t', t") = a(t')a(t") [q A - 2g 2 (q ■ k) 2 + (q ■ k) 4 
U(q,t') = mq,t') + f b 4 (q,t'), 



ft = a(t') 



5q 2 - 2(q • k) 2 



ft 



a(t>) 



3H u\ 



(3.12) 
(3.13) 
(3.14) 



where the inner products are with respect to 8 a b- 
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The ultraviolet divergences of these loop integrals need to be regularized. Previous 
treatments have done this by a process called "adiabatic subtraction" in which one makes 
replacements like |x<j(*)| 2 ~~ * \Xg(t)\ 2 ~ I (*) 1 2 - For sufficiently high q, this difference 
vanishes since the modes are evolving adiabatically, and so all loops become ultraviolet- 
finite. However, this procedure has the disadvantage of affecting answers in the infrared 
regime we arc interested in [7, 8]. Moreover, the intuition is that we are "subtracting off 
the adiabatic behavior," a notion which seems difficult in our problem since the infrared 
modes are in no sense evolving adiabatically (i.e. Xq is n °t even approximately a solution 
of the equations of motion). 

In the appendix we present a simple method for computing the finite effects of the 
loop. If we analytically continue the time arguments to the complex plane, we can take 
the time integrals to run over contours such that along the entire contour, |cj/u; 2 | <C 1 for 
all momenta. This means that we can use the adiabatic mode functions (2.6) for the entire 
time integral. This simplifies the computations and allows us to track which parts of the 
loops come from \ excitations, namely those terms proportional to j3 (recall that a mode 
that is not excited has (3 = for all time). Focusing on these terms also regularizes the 
momentum integrals because j3 ~ e -i 2 / k l % 

G(A ) 



F(Ao) 



0.1 




Figure 3. Dimensionless functions F and G appearing in the answer. F enters into the tadpole 
diagram while G enters into the other diagram. 

The result of these computations is easy to state. The first diagram topology con- 
tributes at most 

l^°M*2^^ F (^)' (3.15) 
while the second diagram topology gives 

wz°v*2^Mr a (Bd- <3 - 16) 

The dimensionless functions F and G appearing here are given explicitly in the appendix, 
eq. (A. 20) and (A. 41); we plot them in fig. 3 for convenience. The key is that they are 

8 This process is actually equivalent to adiabatic subtraction for diagram A, but we note that adiabatic 
subtraction is still only sensible once we have deformed the integration contour so that we actually have 
adiabatic solutions for all time. 
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both O(10 _1 ) and peaked for k ~ a^H, i.e. the effect on the tensor spectrum is dominated 
for tensor modes whose wavelength is roughly the size of the horizon at the time x becomes 
massless. 

In the end, this means that the complete tensor power spectrum, computed to one 
loop and including only the effects from x excitations, is given by (at most) 9 



2 



HM 2 \ \a H J \a H 



(3.17) 



Let us estimate the size of this correction. The energy density contained in the x excitations 
(normalizing a = 1) should be roughly E A = f d 3 q n q u) where n q = |/3| 2 ~ e -i 2 / k i \ s the 
occupation number for the qth mode. We have roughly u ~ q and so this integral gives 
E 4 ~ kf. In order for inflation to proceed, this energy should be subdominant to the 
inflaton potential, so we have E 4 ~ kf <C Vinf ~ H 2 M 2 via the Friedmann equation, thus 
k 2 /M 2 <C H/M. This means that the correction to the power spectrum is bound from 
above by 

k 3 k 

MP « W < 3 ' 18 > 

This quantity had better be quite a bit less than unity since otherwise we are exciting x 
quanta with energy of order the Planck mass. Thus, this correction is quite small. 

4 Classical x backgrounds 

* n 

W // 

* '/ 

• '/ 



II 
II 




II 
II 


n 



II 
II 
II 



_ii 



Figure 4. Some typical diagrams involving background x fields, insertions of which we denote with 
crosses and double lines. Dashed lines are 5\ propagators and solid lines are ( propagators. 

If the quantum fluctuations of a field x become large enough they behave classically. 
We now consider the case whereby quantum fluctuations have been amplified through some 
unspecified process enough to become classical. 

9 The term involving G is precisely the one found in [2]. 
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The slow roll constraints imply that the inflaton potential is much greater than the 
derivative terms of the x held. The corrections to the tensor power from the inflaton are 
proportional to the inflaton potential whereas the corrections from the x background are 
proportional to derivatives of x- The scalar power spectrum depends on derivatives of the 
X background at the same order. Conversely, the scalar spectrum is more sensitive to the 
inflaton potential than the tensors. Thus, the corrections to the tensor-to-scalar ratio r are 
more sensitive to the inflaton potential and are relatively insensitive to the background x 
field. 

The corrections may overcome this suppression given an enormous number of spectator 
fields, however then one must question the use of perturbation theory for studying such a 
system since each field will contribute to the interaction Hamiltonian through at least a 
gravitational coupling. Therefore inference of the inflationary scale H based on the positive 
detection of tensor modes is resilient to the effects of introducing additional scalar fields. 

5 Conclusions 

Detection of tensor modes generated during inflation is currently our best hope of determin- 
ing the energy scale of inflation, which to date is largely unknown. In the simplest models 
of inflation, in which tensor vacuum fluctuations dominate, a measurement of the tensors 
is in one-to-one correspondence with the energy scale. It is conceivable that dynamics in 
more complicated scenarios could destroy this correspondence. 

Corrections to the tensor power spectrum will generally lead to corrections to the scalar 
power spectrum. Recently it has been proposed that sudden excitation out of vacuum 
of some spectator x fields could lead to enhancements of the tensor spectrum without 
significantly affecting the scalar spectrum. 

In this work, we have given a systematic field theoretical treatment of both the excita- 
tion of x an d its effects on the tensor spectrum. Using this framework, we have studied the 
case in which x is a scalar field. We studied the case when the field is excited at a specific 
time, where we have shown that there is a self-consistent upper bound on the corrections 
to the tensor power spectrum. This upper bound constrains the effect to be unobservable. 

These results are consistent with previous studies. It has been suggested recently that 
modifications to this setup, especially the spin of the x field, could change the answers we 
have obtained here. In our language, this is conceivable since it could change the powers of 
the momentum scale fc* appearing in the correction to the tensor spectrum. In scenarios 
with certain vector couplings, gauge quanta can be produced continuously rather than 
suddenly [5]. One could consider the case in which there is a large number of x fields, 
as in axion-based models of inflation [3, 4]; our computation suggests that perturbation 
theory will break down with a large number of these fields, and new methods will have 
to be applied. However, we do conclude that having one (or even quite a few) spectator 
scalar x fields suddenly excited during inflation will not produce a significant effect on the 
tensor power spectrum, and so the detection of tensor modes in these scenarios is still in 
one-to-one correspondence with the energy scale of inflation. 
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Figure 5. Contour used for evaluating /3. 
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A Evaluation of loop integrals 
A.l General setup 

Because the x held evolves non-adiabatically near to, it cannot always be expressed in terms 
of positive and negative frequency mode functions. In computing certain observables this is 
at best a technical headache, but in particular for loop integrals this is problematic because 
it makes it difficult to extract the part of the loop involving excited or "particle" states. 

In order to circumvent these issues, we can extend a trick developed in [11]. The key 
idea is that the non-adiabatic behavior of x is characterized by the fact that |w/o; 2 | 3> 1, 
but only in a small region of time around to- If we analytically continue in t, then we can 
integrate along a contour 7 = j(t) in the complex t-plane such that |w/w 2 | <C 1 for all t. 
Then the adiabatic mode expansion is valid for the entire integral. 

As it turns out, once the integrals are analytically continued in this way, they can be 
evaluated by a steepest descent method. Generally, in evaluating j3 using (2.7) or the time 
integrals in the loop diagrams, we encounter integrals of the form 




(A.l) 
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where / contains some (negative) powers of to and the time-dependence in /, other than 
the frequency terms, is much slower than the time-variation of to. Although to is real and 
positive for all real t, its analytic continuation has two zeroes. Then the integrands contain 
branch cuts from the square root in to as well as (sometimes) poles from factors of to -1 . We 
can thus take the contour 7 along which we perform this integral to be a steepest descent 
contour localized near these points because the exponential is oscillating rapidly except 
very close to the points where to 0. 

Let us set up the analytic continuation and steepest descent contours. Near to and in 
the slow roll approximation we have 



to 



^ + -At-to) 2 -( 3 -^) = k lv(t-t-)(t-t + ). (A.2) 

a a \ z J a 



Here we have assumed that that period of non-adiabaticity is shorter than a Hubble time 
so that a(t) « ao = a(to), and as before k* = At -1 3> aoH is the comoving momentum 
scale of the problem. We have factored to by its zeroes 




t ± = t ±a, * = — (A.3) 



The branch cut(s) of this function must be taken with some care. As described in section 
2, there are different types of behavior as a function of the momentum k. For very low 
momentum modes, a is real and 10 has its zeroes on the real axis. However, for high mo- 
mentum modes, k/aoH > 3/2, a is pure imaginary, and the zeroes of to are symmetrically 
distributed across the real line; these are the modes that affect the tensor spectrum and we 
concentrate on them from here out. When we analytically continue to, we have to choose 
a branch cut to make this a well-defined function, such that to is real and positive along 
the real axis. The contour of integration can be deformed into a steepest descent curve as 
long as we do not cross the branch cuts. 10 

Consider integrals of the form (A.l), with to continued as just described. Let us 
construct a good steepest descent contour 7 along which to evaluate the time integrals. 
First break up the contour as 7 = j a + 70 + jb where 7^ are the contours leading into and 
out of the branch point, i.e. from —00 to t± and from t± to +00, respectively, and 70 runs 
from t± to t near t±. The exponential becomes a product and so 

jj tl (J\ J\ —2* f uxlt" -2i f Ludt" —2* f udt" , . .n. 

at j{Lo(t ),t )e J T" e J ^o e J 'b . (A.4) 



7 



Assuming we can construct a good steepest descent contour, the dt' integral will be totally 
dominated for t? near the branch point t±. Thus the dt' integral will be localized to 70 



10 Many authors write that because the metric's expansion is assumed slow compared to the timescale of 
the excitations, it is consistent to directly set H — in the frequency. This is not quite correct, but the 
phase space of the excited modes is enormously larger than the phase space outside the horizon, and indeed 
the excited modes dominate the loop integrals for exactly this reason. 
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and will get the exponential's contribution from 7 a and 70 only. Moreover, since the time- 
dependence in / is slow compared to oj we can set f(uj(t'),t') = f(uj(t'),to). In other 
words, 

I = I a [ dt'f{u{t')Me~ 2i ^ dt '\ I a = e~ 2l ^ dt ". (A.5) 

Written as such, we have two integrals to do: I a , which gets us from the asymptotic past 
up to the singularity, and the dt' integral near the singularity. To define 7 we need to focus 
on the latter. What we want is a contour such that the exponential is of the form e u with 
u real and negative, so that the integral localizes. Also, we want to take the branch cut for 
uj such that the exponentials will die off transverse to the cut, i.e. we want the cut along 
the curve with \oj\ ~ near t±. It turns out that all of this can be achieved by taking the 
contour straight into the singularity t±, evading the singularity by a small circle, and then 
going straight away from the singularity, see the figure. 

Write t — t± = 5± = r±e l9± , so that we can define the branch of to we want by 
uj = ^/ r _ r+e l ( e -+ e +)/ 2 . Near the branch points, the exponential goes like e u± where 11 

u± = -2i f dt uj » n±5 3 J 2 , n ± = -^4 e±W4 vlM (A.6) 
Jt± 3 a 

The branch cuts should be taken along local angles a± so that the cuts run to ±zoo, so we 
have the restrictions 





(*) 


^**^ 


UHP 


< a + < it 


a + — 2-k<6 + < a + 


LHP 


-ir < a_ < 


a- < 6>_ < a_ + 2-7T 



(A.7) 



on the relevant angles. Now, we want to choose a± such that argu-t is an odd multiple of 
7r/2 along the cut, and such that we can find two rays emanating from S± = to use as the 
steepest descent curve (one ray for the incoming and one for the outgoing direction). It 
turns out this is only possible in the lower half plane. We have &rgu± = &TgQ± + | arg^-t, 
so from the condition (*) we see that in the upper half plane we can only take the cut 
along a+ = tt/2, while on the lower half plane we must take the cut along a_ = —tt/2. 
The steepest descent paths are rays emanating from 5± = along the angles 9± such that 
argu-t = (2n — l)ir for some integers n, ensuring that e u± ~ e~' n± l In the upper half plane 
this means 9+ = (8n — 3)ir/6 and in the lower half plane 9- = (8n — 1)tt/6. However, these 
angles must still obey condition (**) above, so we find that in the UHP we can only use 
n = and thus cannot make a good steepest descent contour. In the LHP, we can use both 
n = —1 and n = 0, corresponding to the rays along 9- = —tt/6 and +7-7r/6, respectively. 
Thus we must take our integration contour in the lower half plane. 

The conclusion of this discussion is that we define the branch cuts with a± = ±tt/2 
(as shown in the figure), and the steepest descent contour is going around the singular 
point t-, consisting of the (incoming) ray arg(5_ = 7ir/6, a small arc around the pole at 
5_ = 0, and the (outgoing) ray arg<5_ = — it/6. With this setup, the integrals near i_ are 

All angles are measured from the horizontal as usual. 
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easy. We expand the factors of oj in / around t = t- and obtain some simple complex line 
integrals. In practice we need the cases / ~ oj/oj = 1/5- (for the f3 integral (2.7)) and 
/ 



oj 



Q~ n 5_ for the loop integrals, where 



Q 



-in/i 



a 



Elementary complex analysis yields the following: 



7<> 



OJ 



h. 



dt'- 



la 



UJ- 



Q 



2m 
Airi 



dt' 



70 



OJ 



QV6(l-i)T ( ) I^T 1 / 3 



(Ai 



(A.9) 



dt'- 



70 



UJ-' 



Q 3 V6(l-i)r( - j |fi-| +1 / 3 . (A.10) 



We also need the integral I a from t = — oo to t = t-. As we move horizontally, oj picks 
up a real part and the exponential just gets a phase. However, going vertically down to t_ , 
the integral gets an imaginary part. Parametrize the vertical coordinate here by s, taken 
from to —a. In our conventions, along the vertical line we have oj = (k*/ao) 2 \/a 2 — s 2 . 
Thus 



exp 



-2i 



i ds ui 



exp 



7T 

2 



3a HY 



(A.ll) 



up to a phase. 

Note that here we have only proven this formula for k S> 3aoH/2. For superhorizon 
modes, the argument goes over almost identically, except that the branch cut is between 
the (now real) zeroes of oj and the contour must be adjusted accordingly. One finds that the 
formula above holds for arbitrary modes, in agreement with mode-matching computations. 

A. 2 Evaluating (3 at late times 

With all of this in place, finding the late-time value of (3 is easy. Starting from (2.7) we 
consider late times t — > oo and take the time integrals along the steepest descent contour 
7 we have just constructed, 



P(k,t-> oo) 



dt 



exp 



2oj{k,t') 

In the notation of the previous section this is just 



-2% 



t' 



dt"u)(k, t") 



(A.12) 



1 



7T 



a-*0 



exp 



7T 

2 



3a HY 
2~k~) 



(A.13) 



Our argument does not fix the phase of /3; this phase is not relevant for many observables, 
in particular the "occupation number" of the x field seen by a late time observer is = 

\m\ 2 - 
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A. 3 Tadpole diagram (A) 

We can extend the technique described in the previous section to compute the loop diagrams 
from section 3. First, for each loop momentum q, we take the time integral to go from 
past infinity t — > — oo to very late times t — > +00 along a contour 7 in the complex t-plane 
such that |^/^ 2 | "C 1 along the entire curve. This ensures that the x mode functions can 
be written in terms of the adiabatic mode functions (2.6). Due to this, we can extract the 
parts of the diagram proportional to (3, i.e. the contribution from x excitations. Defining 
the contour as such, we may deform it continuously to the steepest descent contour 7 
constructed above. After doing the time integrals, we are left with the integrals over loop 
momenta q, which will contain factors of (3 ~ e~ q ' * and will thus converge without any 
need for regularization. For both diagrams our arguments are valid only up to a phase, so 
what we are really going to do is over-estimate all of the contributions by computing the 
modulus squared of each diagram and adding them coherently. 

Let us begin with the tadpole diagram. As there, we label the external graviton's 
(comoving) momentum by k and the internal x propagators by q. Starting from eq. (3.9), 
we take the external time t — > 00, and deform the time integral to the contour 7, so the 
contribution to the correlator is 

A = ~ i j J dt ' h(Q,t')g 4 (k,t')\ Xq (t')\ 2 . (A.14) 

Here f± is given by (3.13), the graviton portion is 

2 H A ., A 2 2k 



g 4 (k,t) = Im [h k (oo)h k (t)j = 2M 4 fc6 [( 1 ~ -j)sinA - A cos A], A = — , (A.15) 
and using (2.5), (2.6) and the Wronskian condition \a\ 2 — \f3\ 2 = 1 for x we have 

\ Xq (t)\ 2 = -L_ f2Re (py/TTfPe-* + 2|/3| 2 + ll . (A.16) 
2crw L V / J 

Let us consider the scaling in momentum (at high momentum) of each of these terms. 
We have ~ q 2 (or slower), co ~ q, and the phase space integral gives q 2 . Thus, if the 
momentum integral is cut off at some scale A, the last term in brackets will diverge as A 

2 / 1,2 

at leading order. However, the other two terms are finite due to (3 ~ e~ q ' *. Moreover, 
we are trying to study the effect of the presence of excitations, i.e. non-zero f3. Therefore, 
we will discard the first term from here out and focus on the finite contributions from the 
latter two terms. 

Consider the integral from the first term in (A.16), 12 

A l = ~ i 7W% / ^ m,t')g,{k,t l )P{q,t'Ul+ (3(q,t') 2 . (A.17) 

Here /3 = f3(q,t) = for t < to an d becomes non-zero at to- Again assuming that the 
timescale of non-adiabaticity is much shorter than the Hubble scale, we can take the time 



2 We first focus on the terms from j% because the exposition is simpler; the results from f\ are identical. 
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arguments of /J 1 , 54, and a to be to since they are all varying much more slowly than 00, 
and so we just need the time integral I\ from the first section of the appendix. Explicitly 
(dropping subleading terms in /3 <C 1) 

7^3/4(9, to) 9i (q, toMq, to) x ^ (A.18) 
where Aq = 2k/aoH and F is the dimensionless O(10 _1 ) function 



^(Ao) 



1 



A 2 

1 r ) sin Aq — Aq cos Aq 



(A.20) 



Note that this function is peaked for k/a$H ~ 1, i.e. for graviton modes fc exiting the 
horizon at t = to, see fig. (3). 

Now we need to do the loop integral involving the second term in (A. 16). This term 
does not have an oscillating phase but can be easily bound from above. To define the 
integral in the first place, we do the same trick as before, taking the time integral to run 
along a large semicircle in the lower half plane so that we can use the adiabatic mode 
functions. Having done so, consider the integrand 

A / , ^/l(')9l(*),o,_M2 



This integrand has no branch cuts, but still has the poles where w = 0. Now, we may deform 
the contour onto the real time axis, plus a circle around this pole. The contribution from 
the pole vanishes as we shrink the circle, since near the pole uj ~ 5 l l 2 so the integrand 
scales with the circle's radius as ~ R 1 ' 2 -> 0. Thus we only have to do the integral on the 
real axis. Again (3 is zero until t > to. After to, the integrand is bounded from above by 
its value at to (uj is increasing and the other functions are redshifting) , and decays to zero 
within a few multiples of the non-adiabaticity timescale At, so we can overestimate the 
time integral by taking the value at to and multiplying by At = A;" 1 . We can then perform 
the momentum integral. The result is 

WSOWx^fW (A.22) 
exactly the same as the other term in the diagram, up to a numerical coefficient. 
A. 4 Diagram B 

Finally, we need to evaluate the second diagram topology. There are four diagrams here: 
RR, LL, RL and LR. We will use the same techniques from the previous two sections, 
although this diagram introduces the additional complication that it has two internal time 
coordinates t' and t" . However, the loop integrand is peaked along the diagonal t' ~ t" 
due to redshifting and relative phases. Once we reduce the integrals to ones over only t' , 
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we can apply the steepest descent method to localize them to t' ~ to as in the previous 
section. 

With this in mind, we start with eq. (3.10) and (3.11). Once again taking the external 
time t — > oo and taking both time integrals along a large contour in the complex t',t"- 
planes, we see that the contribution from the RR diagram plus the LL diagram gives 

B RR+LL = J dt '^ Q Mq, t', t")g 3 (k, t', t")h 3 (q, q', t> , t") (A.23) 

where here f 3 is given by (3.12), the momenta satisfy k + q + q' = 0, g 3 comes from the 
graviton propagators 

g 3 (k, t', t") = Re (Moo) 2 W)W)) (A.24) 
and h 3 comes from the \ propagators, 



h 3 (q,q',t',t") = X g (t')x q (t")x q '(t') Xq '(t"). (A.25) 

The momentum integral is going to be dominated by the internal momenta q ~ k* while 
the external momentum is only relevant around k ~ a^H <C fc*. Thus we may approximate 
q' ~ q. With this approximation, let us insert the adiabatic mode functions (2.5), (2.6) 
into (A.25). We obtain 

h 3 ( q ,t',t") = 4^)^)^,)^) [ 4 I«| 2 |^| 2 

+ 2\a\ 2 a(3e- 2ie ' + 2a/3|/3| 2 e+ 2ie ' + 2a/3|/3| 2 e - 2ie " + 2|a| 2 a^e +2ie " 
+ | a | 4 e- 2iAe + |/3| 4 e +2iAe 
+ a 2 (3 2 e- 2ie 'e- 2iAe + a 2 f e +2ie ' e +2iAe 

Here the phases are 

/t' ft" ft" 

cLtuj{t), 6"= / dr u(t), A9 = 6' - 9" = / dr w(r). (A.26) 
-oo J —oo J t' 

This somewhat fearsome expression for h 3 shows that there are four types of terms, 
organized by line. The first term has no oscillating phase but is straightforward to estimate 
using the same logic as A2 above. It gives a contribution 

! ^ f dt'dt"d 3 qf 3 (q,t',t")93(k,t',t") ] . 2]R]2 ( , 

Brr+ll- I (27r)3 a^t'W^it'Mt") 1 1 Wl ' { ' 



As explained earlier, the t" integral will localize to t" ~ tf, so we may set t" = t', J dt" — > At. 
The remaining integral over t' has zero integrand for t' < to an d the integrand is bound 
from above by its value at t! = to (again because the frequency increases and the rest 
of the functions are redshifting) . Thus we can set all the time arguments to to and take 
J dt' — > At, obtaining an easy momentum integral which yields 

\B RR+LL \ < 0(1) x 2 M2 k3 h % GrR + ll(\o) (A.28) 
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which has the same parametric dependence as the previous diagram, but a different dimen- 
sionless function from the graviton propagators, 



A o 



A 2 

1 — ] cos Ao + Aq sin Aq 



(A.29) 



Here once again Ao = Ikja^H. Although this function is divergent as Ao — > 0, we will see 
that the total contribution from all four diagrams gives a finite result. 

The second type of term goes over similarly to the previous diagram. Consider for 
example the first such term, which has the phase ~ exp(— 2i f lo), just like the oscillating 
terms from the previous diagram. As argued earlier, the t" integral is dominated near 
t" ps t' so again we set t" = t',f dt" — > At. Other than the phase 9' = J** co all other time- 
dependent functions are varying slowly, and we can localize the t' integral to a steepest 
descent contour near t' = t- as before, so we get (again dropping subleading terms in j3) 

/dt'd?a e~ 2i f u 

This time integral is I2 from above. Using that result we are left with an easy momentum 
integral that yields 

\B\ R+LL \ < 0(1) x 2 M2 k z hM* GrR+Ll{X ^ (A - 31) 

just like the first term. 

Terms from the third line require a bit of extra work, although they actually give the 
same result. We need to be a little delicate with the t" integral. This should still be 
dominated by t" ps t' , but now we note that in this approximation 



/" 

j —00 



dt „ e -2iAe = [ t dt „ e -2itf dru{T) w , A 32) 

Loo Mt') 



After estimating the t" integrals in this fashion, we have reduced everything to the case 
from the previous paragraph and again we obtain answers of the form 



\B RR+LL \ <0(1) x 2 M2 k 3H k L Gl!, "- LL{X{l) - (A " :W) 



Finally, to evaluate the terms in the last line, we essentially combine the arguments from 
the previous two types of term. Namely, the t" integral is replaced by J dt" — > i/2u;(t'), 
leaving us with a time integral of the form 



/ 



e -2i J uj 

dt'-^y = I 3 . (A.34) 



Using this and the usual arguments about t' « to we again obtain momentum integrals 
which evaluate to 

H 2 A- 3 

\B RR+LL \ < 0(1) x — —-^ Grr+ll (\ ). (A.35) 
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In summary, what we have found so far is that all the terms from the right-right and 
left-left diagrams give the same type of contribution, 

\B RR+LL \ < 0(1) x 2 ^3^ 2 Gwu(Ao). (A.36) 

Now, we need to do the right-left and left-right diagrams. These are in fact identical so we 
only need to do one. Writing out the integrand starting with (3.11) and following the same 
steps used to simplify it as in the previous case, one sees immediately that the integral is 
exactly the same as the RR case, except that the graviton propagators are different 



93 -> h = \h k (oo)\ 2 h k (t')h k (t"). (A.37) 
The effect of this is that the combined diagrams give 

\B RL+LR \ < 0(1) x 2 * 2 2k3 J^ Grl+lr (\ ), (A.38) 

where 

G RL+LR {\ ) = ^{l + (A.39) 

and the coefficient we have dropped is the same as the RR + LL case. The relative negative 
sign between the RR and RL diagrams (3.10), (3.11) then conspires to make the total result 
finite at Ao — > 0, and we find that the total contribution to the tensor two-point function 
from the B diagrams is 

WOll)"^^ (A.40) 
where again Ao = 2k/aoH and G is 

12 



G(A ) = T3 



A 2 / A 2 \ 
1 + ^ - f 1 - y J cos A - A sin A 



(A.41) 



This function is 0(10 1 ) and peaked for k ~ a^H/A, similar to the function -F(Ao) appearing 
in the other diagram; see figure (3). 
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